Smooth quantum— classical transition in photon subtraction and addition processes 
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Recently Parigi et al. [Science 317, 1890 (2007)] implemented experimentally the photon subtrac- 
tion and addition processes from/to a light field in a conditional way, when the required operations 
were produced successfully only upon the positive outcome of a separate measurement. It was veri- 
fied that for a low intensity beam (quantum regime) the bosonic annihilation operator a does indeed 
describe a single photon subtraction, while the creation operator describes a photon addition. 
Nonetheless, the exact formal expressions for these operations do not always reduce to these simple 
identifications, and in this connection here we deduce the general superoperators for multiple pho- 
tons subtraction and addition processes and analyze the statistics of the resulting states for classical 
field states having an arbitrary intensity. We obtain closed analytical expressions and verify that for 
classical fields with high intensity (classical regime) the operators that describe photon subtraction 
and addition processes deviate significantly from simply a and a* 1 . Complementarily, we analyze in 
details such a smooth quantum-classical transition as function of beam intensity for both processes. 

PACS numbers: 03.65.Ta; 42.50.Ar; 42.50.Lc 



I. INTRODUCTION 

It is well known [l| that the probability for absorbing 
one photon per unit time from an electromagnetic field 
is proportional to the average value of the ordered prod- 
uct of the negative and positive frequency electric field 
operators over some quantum state p, also known as the 
statistical operator. In the simplest case of a single-mode, 
that probability can be expressed in terms of the stan- 
dard bosonic 'annihilation' and 'creation' operators a and 
at, satisfying the commutation relation [a, at] = 1, as 



P = 7T1- (apa 



(1) 



where p stands for the field state just before absorption 
and 7 is an appropriate coefficient. After interacting with 
some detector that absorbs one photon, the field makes a 
transition to a new state, which can be formally described 
by the action of a photon subtraction superoperator (PSS) 
V as 



P 



= p- l v P , 



(2) 



where p' represents the field state immediately after the 
subtraction of one photon and P = Tr(Dp) is the proba- 
bility for that process. 

The hcrmiticity of p' is always assured whenever T> 
takes the form 



Vp = jOpO^ 



(3) 



where O is some operator responsible for the subtraction 
of one photon from the field and the explicit form of O 
depends on the details of the field-detector interaction. 
Since the 1960s several models were proposed d, 0, H, @ 
(see Q for more references therein) and the first one [2| 
used the rather simple identification O — a, when Eq. 
© becomes 



and we shall refer to it as A-model. Although such a form, 
Eq. ([5]), seems quite natural in view of equation fT]). this 
choice was, as a matter of fact, intuitive, although later, 
under certain assumptions - such as the weak coupling, 
low intensity and short interaction time - the authors of 
Refs. [5, 8] were able to derive it from doing a microscopic 
analysis. 

Nonetheless, if those assumptions are replaced by oth- 
ers, one can obtain different superoperators T>. A family 
of PSS based on the nonlinear lowering operators of the 
form O = (1 + n) _/3 a, where h = a^a, was derived in Ref. 
[8j. The special case (3 = 1/2, called E-model, 



£p = E_pE + , E- = (l + n) 



-1/2, 



(5) 



Ap = japat 



(4) 



was originally proposed ad hoc in (9l.[Tfj|. Later, analy zing 
carefully a microscopic model, it was shown @, M, [If 
that both superoperators A and £ are particular cases of 
a more general one, where A {£) is specific for a small 
(large) mean photon number. 

Recently, the A-model was subdued to an experimental 
verification in the low photon number regime and weak 
field-detector coupling [1 31 ] . That experiment employed a 
beam-splitter and a single photon on/off detector (SPD), 
within a simple scheme as illustrated in the Fig. [T^: a 
low intensity beam, with a small mean photon number 
per unit time, is prepared to hit a low reflectivity beam- 
splitter. The reflected beam is continuously measured by 
the SPD, and whenever it clicks, the transmitted field 
statistical operator can be described, approximately, by 
the A-model 0. 

However, as already observed in [11] , A is an un- 
bounded superoperator and some physical inconsistencies 
appear [l4[. For example, A does not hold for large values 
of the initial mean photon number, since the probability 
of a photon subtraction, P = jTr(apa^), may become 
larger than 1. Moreover, applying A on some specific 
field states, the calculations predict results which look 
counterintuitive. For example, operating with A on the 
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FIG. 1: (a) Scheme for conditional photon subtraction by 
means of a beam-splitter (BS) and a photodetector (D). (b) 
Scheme for conditional photon addition using a pump beam, 
a non-linear optical crystal and a photodetector. 

thermal state, which has the photon number distribution 



where no is the mean photon number, for post-selected 
state we obtain (n)^ = 2hq. Although this prediction 
was confirmed experimentally in the 'quantum limit' of 
small mean photon number, no < 1, it is counterintuitive 
from the classical point of view, namely, a high intensity 
field cannot double the intensity due to a one-photon de- 
tection. So A can describe satisfactorily the low intensity 
field but not a high intensity one. Quite differently, us- 
ing the E- model we get (n)g = no, which is sound for a 
high intensity field after one-photon detection, whereas it 
does not reproduce the observed doubling for low inten- 
sity. Therefore, it is interesting to analyze the behavior of 
the PSS in the 'classical limit' (large mean photon num- 
ber) and how the pattern of post-selected state statistics 
changes smoothly and continuously as function of the 
field intensity. The quantum- classical transition in pho- 
ton subtraction process regards the passage between the 
two extreme situations for classical field states that can 
be prepared with arbitrary low and high intensities. 

Besides the photon subtraction, one may also imple- 
ment the photon addition operation by the conditional 
stimulated down-conversion in a non-linear optical crys- 
tal, as reported in [l3| and illustrated in the Fig. QJ>. 
Inside the crystal the pumped photons may decay spon- 
taneously into two entangled photons with energies that 
sum up to that of the parent one. Upon detecting one 
of these photons (known as the trigger photon) along a 
particular direction, the other photon state is unambigu- 
ously determined. If one injects a seed light into the 
crystal, the stimulated emission may occur into the same 
mode, and the detection of a single trigger photon (by 
means of a SPD) indicates the conditional generation of 



the photon-added state. Using this scheme, the one pho- 
ton addition superoperator (PAS) for low beam intensity 
is 

Ap = -f'a^pa (7) 

(we call it the A + -model), which was experimentally im- 
plemented for a low mean photon number [l3| . Since the 
superoperator A is unbounded, with the probability of 
photon addition becoming larger than 1 for large field 
intensities, in analogy to the E-model we also define the 
E + -model as 

:.p /■:./</•: . (8) 

that will also be analyzed and discussed below. 

In the present paper we generalize the theoretical anal- 
ysis presented in [1 31 ] for both, the photon subtraction 
and photon addition operations, by substituting the SPD 
by the fc-photon resolving [l9j or nonresolving [20| de- 
tector, and obtain formally PSS and PAS valid for an 
arbitrary field at any intensity. Although the resulting 
formal expressions for the PSS and PAS are compact, in 
general, it is not straightforward to obtain closed analyti- 
cal expressions for these operations for an arbitrary initial 
field state. Therefore, it is desirable to have some sim- 
ple approximate expressions for the PSS and PAS that 
hold in a specific regime of parameters and can be easily 
evaluated for an arbitrary field state, providing a sim- 
ple means of predicting the outcome of the experiment. 
In this connection, we show that in the quantum limit 
(small photon number), the A- model and A + -model are 
good approximations to the exact PSS and PAS for any 
field state, respectively, while in the classical limit (large 
photon number), the E- and E + - models are more appro- 
priate for a nonresolving detector and the 'classical' field 
states: coherent, thermal and 'mixed light'. 

The paper contains three additional sections. In sec- 
tion[TT]we obtain exact analytical expressions for the PSS 
and analyze the behavior of the photon subtraction prob- 
abilities and post-selected states statistics for any inten- 
sity of the field and for different kinds of 'classical' field 
states: (a) coherent, (b) thermal and (c) 'mixed light' 
[2l| . which, in principle, can be produced in the labo- 
ratory with arbitrarily low or high beam intensities. In 
section UlTl we do the same for the PAS, considering (a) 
coherent and (b) thermal field states. In section IIVI we 
present a summary and our conclusions. 

II. PHOTON SUBTRACTION 

After passing through the beam-splitter a small frac- 
tion of the incident field (signal) is reflected into the mode 
b, initially in the vacuum state |0(,)(0{,| (Fig. QJi). The 
post-selection procedure - in which k > photons are 
detected by a detector placed in the reflected path - al- 
lows to express the transmitted field state as Eq. 
where the exact PSS is written as 

V k p = Tr b [M k U(p®\O b )(O b \)Ui] , (9) 
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where p is the incident field state. Here 
U = cxp [9 (a)b - atf)] 



(10) 



is the beam-splitter operator, where a (signal beam) and 
b (reflected beam) represent the two modes resulting from 
the incident beam, and 9 is a parameter related to the 
reflectivity (R = sin 2 9) and transmittivity (T = cos 2 9) 
coefficients. The operator M k — Yl'iZk ^ il'&X^&l) that 
acts on the reflected mode, stands for the action of the 
detector (referred as T>k, for short). Tj — 1 for a nonre- 
solving fc-photon detector (NVk), that clicks whenever 
k or more photons are absorbed [H, and T; = 5i ik 
for the resolving fc-photon detector (RDk), that clicks 
when exactly k photons are absorbed [19( . The SPD cor- 
responds to NT>1 according to this notation. 

The resulting PSS is a generalization of the Eq. ([3]) 
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v kP = J2r l (i b \u\o b ) P (i b \u\o b ] 



t 



(11) 



l=k 



The operator [7 has the form exp [—9 (K + + K-)] where 
K + = b^a, K_ = —ba^; a third operator, K = (b^b — 
a^a)/2, is necessary to close the su(l,l) algebra, they 
constitute the generators of the SU(1,1) group. Using 
the factorization theorem [IH, [l(| [l?], [HI we can write 
OH as 



l=k 



(tan 2 9) 



which is the exact and complete PSS compatible with 
the experimental setup reported in For a small 

mean photon number no and R <C 1 one has roughly 



Tr 



(cos 9) n p (cos9) n 



obtains 



(1 — noR), so for uqR <C 1 one 



V k p~AkP=—" k »^ k 



R k 

-M a " pa ' 



(13) 



for both kinds of detectors. For k = 1 one retrieves the 
A-model, whereas for A: > 1 we have a straightforward 
generalization of A-model for multiphoton subtraction. 

The photon number distribution of the post-selected 
state for the k photons detector is 



l>n,Vk ~ r Vk "n i 



where 



ef fc) =E T ' 



l=k 



T n R' Pn+ i 



(14) 



(15) 



and P vk = 
subtraction, with p 
distribution. 



50 Q 1 ^^ is the probability for fc-photon 



denoting the initial photon number 



Below we specialize to three different states of the field, 
which are currently produced in the lab: the coherent, 



FIG. 2: (Color online) NVl and RV1 results for the 'mixed 
light' state with mean photon number no compared to the 
predictions given by superoperators A and £. Here n c = n t /4. 



thermal and the 'mixed light' [21| . The latter is a 'mix- 
ture of the thermal and coherent radiation' (2^|, whose 
the photon number distribution is 



■a 1 (cos Of p (cos Of a) 1 . (12) (mi) 



exp 



l + n t J (n t + 1) 



n+l 



Li, 



nt ti- 



nt) 
(16) 

where n c and n t are the mean photon numbers of the co- 
herent and thermal part of the mixed light, respectively, 
uq = n c + n t is the total photon number and £„(•) is a 
Laguerre polynomial. For these three states the resulting 
expressions for p' nT>k , Pr> k and the two lower moments 
of the photon number distribution for the post-selected 
state are given in the appendix IA1 

For a small mean photon number, hqR <C 1, the ex- 
pressions for the photon subtraction probability and the 
two lower moments of photon number in the resulting 
state [see Eqs. in appendix [X] are approximately equal 
to the expressions corresponding to the generalized A- 
model. Thus, in the quantum regime (low intensity held), 
A k is a good approximation for the PSS, and the detec- 
tion of one or more photons may significantly increase 
the mean photon number of the post-selected state. On 
the other hand, for a high intensity field, n Q R 3> 1, the 
results are quite different from those predicted by using 
A k . This behavior can be appreciated looking Figs. [5] 
and[3l where we have plot Pdi, { n )vi an( l ( n ( n — l))x>i> 
together with the corresponding predictions of A\ and 
£\, for the mixed light state with n c = n t /A (Fig. \2§ and 
n c = I0n t (Fig. [3]), setting R = 10~ 2 . One can perceive, 
from the figures, that for n^R ^> 1 the probabilities and 
the factorial moments for NVl are close to those pre- 
dicted by the generalized E-model. Therefore, for the 
considered classical field states, in the classical regime £ k 
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FIG. 3: (Color online) Same as Fig. [2] for n c = 10n t . 

represents better the photon subtraction process (using 
a nonresolving detector) and, as verified empirically, the 
mean photon number does not increase upon the pho- 
todetection. The smooth transition between the A- and 
E- models occurs in the region uqR ~ 1. 

The important effect of the measurement back-action 
may be clearly verified by comparing the outcomes of 
the instantaneous fc-photon detection (using either NT>k 
or RDk) to the sequential detection of fc photons (one 
by one, as a sequence of discernible clicks), using an 
array of k SPD's and detecting one click in each one. 
For the sequential k photons detection the PSS is S k = 
(XVi) fc ,where Pjvi is given by Eq. (fT2")l with T; = 1. 
The resulting expressions are given in the appendix [B] 
and plotted in the Fig. |4j Looking at it we can compare 
the outcomes of the three kinds of detectors, S2, NT>2 
and RD2, for the thermal state and R = 10~ 2 . It turns 
out that for a sequential counting, the probability of de- 
tecting k photons is always higher than that calculated 
by admitting an instantaneous detection. Moreover, for 
rio-R "C 1 the mean photon number in the post-selected 
state is always higher than for instantaneous detection, 
while for n$R > 1 different mean values are predicted, 
such that (n) Nvk > (n) Sk > {n} RVk , see Fig. @] 



III. PHOTON ADDITION 

The fc-photon addition superoperator is defined as [l3j 
B k p = Tr b [M k up\0 b )(0 b \vl] , (17) 
where M k = J2"iZ k ^i\h)(lb\ ( as defined previously) and 
u = exp [\(ab-aW)] (18) 



FIG. 4: (Color online) Comparison of 2 sequential detections 
with SPD's to the 2-photon detection using either the resolv- 
ing or the nonresolving detectors, for the thermal state. 



is the operator describing the parametric down- 
conversion process with gain factor A <C 1. The op- 
erator u contains two generators of the su(l, 1) alge- 
bra, K+ — a^tf, K- = ab, while the third one is 
K = {a^a + tfb + 1) /2. So Eq. (JT7J) becomes [3 

B kP = tf> ( tan ^ A )' Q tz ( cosh -i X f p ( cosh -i X f 

l=k 

(19) 

and we define r = sinh A and t = cosh A. For a small 
mean photon number one obtains 

D feP ~ A kP = r -^a) k pa\ (20) 
k\ 

which is the generalization of the A+ -model. As like 
as A k , the superoperator is unbounded, with the 
probability of photon addition becoming larger than 1 
for no 3> 1. Therefore, the A + -model cannot stand for 
high field intensities. Here we shall study the behavior 
of the lowest photon number moments in the transition 
from the quantum regime, hq -C 1, to the classical one, 
rio ^> 1, and compare the results with the predictions of 
the generalized E + -model, M k p = E\pE^. Contrarily to 
Afe, superoperator is bounded, presenting probability 
P&k = 1 and the mean value (n)„, = hq + k for fc-photon 
addition operation. 

The photon number distribution of the fc-photon added 
state is given by 

„' _ p-lo(Dfe) foil 
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FIG. 5: (Color online) iVBl and RBI results for the thermal 
state with mean photon number no compared to the predic- 
tions of superoperators A and E. 



where 

er)=^T / ^Ji"+Vp„_ ; (22) 

and the probability of the fc-photon addition is Pjjfc = 
Y^^Lo 0n Dfe ' i independently of the number of photons 
present in the field. In order to make the physics more 
transparent, we assume the simplest situation k = 1, for 
which P A1 = r (n + 1) and (n) m /n a = 1 + Uq . 

For the two different field states, the coherent and the 
thermal, we get the expressions given in the appendix [Cl 
For ngr <C 1, in both cases, the expressions become very 
close to those predicted by the A + -model, so in the quan- 
tum regime the photon creation operator describes 
accurately the photon addition process. In the regime 
uqt 3> 1, for the coherent state the predictions of the 
A+- and E+- models for the moments (n) and (n (n — 1)) 
are very similar, while for the thermal state, according to 
the A_|_-model, the mean photon number is roughly twice 
the one predicted by the E + -model. In Fig. Owe show 
the behaviors of the photon addition probability P, (n) 
and (n (n — 1)) for the thermal state as function of nor 
for A = 1CP 2 . We see that by using a nonresolving de- 
tector for ngr 3> 1, the expressions approach the results 
of the E + -model, so in the classical regime the A + -model 
ceases to represent the PAS and the E + -model becomes 
more appropriate. The transition between the A+- and 
E+- models occurs in the region nor ~ 1. Finally, notice 
that, although the photon resolving and nonresolving de- 
tectors are quite similar in the quantum regime, they are 
completely different in the classical one, as clearly seen 
from Figs. H El Hand El 



IV. SUMMARY AND CONCLUSIONS 



In summary, we have analyzed the quantum-classical 
transitions for the photon subtraction and for the pho- 
ton addition processes from or into the field. We used 
the usual formal representations for the beam-splitter, 
the optical nonlinear crystal and the /c-photon resolving 
or nonresolving detector, as it was described in [l3| . We 
considered three classical field states - coherent, ther- 
mal and the 'mixed light' - to illustrate the transition 
from quantum to classical regimes and studied the pho- 
ton number factorial moments in the post-selected state 
as function of the intensity of the prepared field, that goes 
on either a beam-splitter or a nonlinear crystal. We ob- 
tained the formal expressions for the PSS and the PAS, 
valid for any input state, and for the considered clas- 
sical states we derived closed analytical expressions for 
the photon number distribution of the post-selected state 
and the associated lower factorial moments, as well as the 
photon addition and subtraction probabilities. 

We found that in the quantum regime (small pho- 
ton number) the PSS can be described approximately 
by the generalized A-model, while the PAS by the gen- 
eralized A + -model for any field state. The mean pho- 
ton number in the post-selected state may increase sig- 
nificantly due to the photon subtraction. On the other 
hand, in the classical regime (large photon number) the 
generalized A- and A+- models lose the validity and, for 
a nonresolving photodetector and the considered classi- 
cal states, the PSS (PAS) is better approximated by the 
generalized E- model (E + -model), and the mean photon 
number necessarily decreases upon a photon subtraction. 
Thus, one may associate P_ used in to the classical 
photodetection operator, in the same way as the bosonic 
annihilation operator a is associated to the quantum pho- 
todetection operator. 

In conclusion, the exact expressions of photon subtrac- 
tion and addition superoperators, Eqs. (fP2|) and (|19p . 
cannot always be written as jOpO' or jO^ pO, they re- 
duce to these simple expressions only in some limits of 
the field intensity or for a photon-number resolving de- 
tector. More importantly, they depend on the form how 
detection is done, so the post-selected field state will de- 
pend essentially on the way the experimenter chooses to 
probe it. 
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APPENDIX A: EXPRESSIONS FOR PHOTON 
SUBTRACTION 

(a) The coherent state |a) is insensitive to the detector 
outcome, since it is an eigenstate of the lowering operator 
a appearing in Eq. (fT0|) . so the resulting post-selected 
state is still a coherent state |acos6>), although having 
the field intensity attenuated. We get for the photon 
number distribution of the post-selected state 

Pn,NVk = P'n,RVk = ^ W f /«!. (Al) 

The normalized factorial moments in the post-selected 
state are the same for both detectors, 



(n) T>k /n = T, {n{n-l)) vk /nl = T 



(A2) 



although the subtraction probabilities associated to the 
detectors are different 



n - e 



_ x n + Rn\ 
(1 + Rn t f 



(A10) 



(n (n - l)) NVl = P NV1 T 2 {n 2 , + n t (n + n c ) (All) 

4- 



_ r nl + 2n c n t + n? (l + 4Rn + 2R 2 n 2 t ) 



(1 + Rn t f J 

where x = Rn c j (1 + Rn t ). For the RVk the correspond- 
ing expressions are quite simple for any fc 

P R v k = e- x - {R " t)h k+1 L k (A12) 



(1 + Rnt) 



iRVk 1 + Rnt 



Lk-i 



(A13) 



t^RVk — e u — — — , -nvpjt 



1-e 



-UqR 



^ (n R) 1 



^ l\ 

1=0 



(A3) 

(b) For the thermal state the photon number distribu- 
tions of the post-selected state for the k photons resolv- 
ing/nonresolving detector are 



Pn,RT>k 



Pn,NVk 



n + k\ ( 1 + n R 
n J \ 1 + n 



fc+l 



n T 
l + n 



1 + n R 

HqR 



(noTf 



k-l 

£ 



n + I 



n 



(1 + n T) 
n R 



71+1 



(A4) 



(A5) 



(n TY 



,_„ v ■ - vl + ?V (1 + no) 
The fc-photon subtraction probability is 



n+l 



(noii)" 



(l + n R) K 

and the normalized lowest factorial moments are 
( n )vk + H vn R 



T- 



l + n R 



(A6) 



(A7) 



(n (n-l)) vk _ t2 (1 + fc) (2 + k) + 2im fl (2 + fc + np-R) 
»o ~ ' (1 + /'-,/.' 

(A8) 

where v = 1 for the TVDfc and i> = for the RVk. 

(c) For the 'mixed light' and the nonresolving detector 
NT>k, the resulting expressions are quite extensive, so we 
write out explicitly only the expressions for k = 1, 



Pnvi — 1 



1 + Rn t 



(A9) 



(n (n — 1)) 



RVk 



n t T 



{z(l+z)-y(A14) 



1 + fln 4 
2kzL k ^i fc (fc - 1) L fe _ 2 



where Lfc = Lk(y), y = -n c / [n t (1 + Rn t )} and z 
l + 2k-y. 



APPENDIX B: EXPRESSIONS FOR 
SEQUENTIAL PHOTON SUBTRACTION 

For the sequential counting, in the simplest case, fc = 2, 
we obtain: 

(a) For a coherent state p' n S2 = e~ n " T (n T 2 ) /nl and 
P S 2 = e - nnR ^ +T '> (e noRT - 1) (e noR - 1) 

( n )s2 _ T 2 ( n ( n - 1 ))S2 _ rj^i 



rpl 



n 



(Bl) 



(b) For a thermal state 

(n T 2 ) 



Pn,S2 = 9l 



(n T 2 Y 



(l + n T 2 ) n+1 (l + 7i T)" 4 
(n T 2 )" (n T 2 ) 



T (B2) 



[l + n (l-RT)] n+1 ' (l + n ) 



\n+l 



and 

Ps2 = 9\ , 
where 



{n) S 2 ^ T 2 92 {n(n-l)) S2 = ^ 



Hi) 



91 



Hi 



9i 
(B3) 



g n = + n R] n -[l + n RT] n +[l + n R (1 + T)} r 

(B4) 

These expressions are quite different from the corre- 
sponding expressions obtained in the appendix [Al 
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APPENDIX C: EXPRESSIONS FOR PHOTON 
ADDITION 



-P/VBl — 



rt(l + n ) 
1 + n^rt 



(C9) 



(a) For a coherent state and the one-photon nonresolv- 
ing detector we get 



Pn,NBl — ATBl te 



( ^» T ( U °\ ("0*)" 



Pnbi =l—te 



(a) 



A. L 



1 + r + r/n - i 2 e~ rt ™° 



n 



Pi 



Ml 



(CI) 



(C2) 



(C3) 



(n (n-l)) Nm = t- 2 + 4r/(n t) + 2r 2 /n 2 - t 3 e~ rt "» 

where L n {-) is a Laguerre polynomial. For the A+- and 
E+- models the resulting expressions are 



(n) A1 /n = 1 + n Q 1 + (n + 1) 
(n(n-l)) A1 /n 2 = 1 + 471^ 
(n (n - 1)) E1 /n 2 , = 1 + 2n~ 1 



(C5) 
(C6) 
(C7) 



(b) For a thermal state and the one-photon nonresolv- 
ing detector one has 



Pn,N01 — -P/VB1 * 



(rc *) n 



(rt + n )" 

. \n+l /-, . \7l+l 

(l + n ) (l + "o) 



(C8) 



(«>atbi _ i^ 1 + r /"o) - [*/ (nort + I)} 2 



n 



Pi 



NI:1 



(CIO) 



(n(n - 1)) WD1 _ „ (t- 1 + r/n ) 2 - [t/ (n rt + l)f 



For the one-photon resolving detector we get 

(l + n Q rt) 2 {n t) n 



(Cll) 



Pn,HBl 

Pr»i 

( n )_RDl 

(n(n - l)) mi 



= 2t 



tn (l + n ) 
rt 2 (1 + np) 
(l + n rt) 2 
1 + t + rtp 1 

1 + nort 

2 + t + 2n 1 
(I + n rt) 2 



— (C12) 
(C13) 
(C14) 
(C15) 



The expressions for the A+- and E+- models are 

(n) kl /no = 2 + no" 1 (C16) 

(n(n-l)) A1 /n 2 - 6 + 4n„ 1 (C17) 

(n(n-l)) E1 /r^ = 2 + 2n 1 . (C18) 
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